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ABSTRACT 

We present a new approach for calculating liquid 
water content in growing atmospheric ice (LWCi) from in- 
situ measurements of both capacitance and conductance. 
Computer simulations and experimental data analysis of 
the four different connections between liquid and solid 
parts of atmospheric ice show that a combination of series 
and parallel connections is the most appropriate model for 
calculating LWCi. Assuming this combination model, and 
using in-situ measurements of both capacitance and 
conductance, we were able to accurately determine the 
time-variation of LWCi. 

1. INTRODUCTION 

The physical properties of atmospheric ice play a 
crucial role in the icing of structures such as aircraft, 
helicopters, power lines, and towers (Minsk, 1983). Since 
the presence of unfrozen water significantly affects these 
physical properties, the liquid water content in growing 
atmospheric ice (LWCi) is one of the most important 
parameters to measure. While intensive theoretical and 
experimental studies of atmospheric ice have been done 
(e.g., Bragg et al 1997, Druez et al 1986, Lozowski et al 
2001), most were conducted after ice formation was 
complete and therefore do not provide data about physical 
properties during the dynamic phase change period. Thus, 
an in-situ study of LWCi is necessary to fully understand 
the physical properties of atmospheric ice. 

Our approach for measuring LWCi makes use of the 
well-known fact that ice and water differ greatly in their 
electrical properties such as their dielectric constant and 
their electric conductivity (Petrenko and Whitworth 1999). 
Thus, if s possible to distinguish between water and ice by 
measuring these electrical properties. This paper describes 
our approach in detail for using in-situ measurements of 
both capacitance and conductance to calculate LWCi. 

2. EXPERIMENTS 

An icing wind tunnel was used to simulate 
atmospheric icing (Fig. 1). Ice was grown with a fine 
water-spray under turbulent-flow conditions in a 20 cm x 


20 cm test section of the wind tunnel. Water droplets of 
distilled water were produced by a spray nozzle, 
manufactured by Spray Systems Co. having air cap 
#173160 and fluid cap # J2050, at an air pressure of 140 
kPa and a water pressure of200-350 kPa. 



Figure 1. Schematics of experimental setup. 1: Di-water 
tank (300 ml), 2: HP4192A, 3: PC, 4: Air compressor. 



Figure 2. The copper-foil grid electrodes on the ice 
collector (Plexiglas). The electrodes had a width of 0.68 
mm and were spaced at 0.59 mm. The electrode thickness 
was 0.075 mm. 
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To accurately determine LWCi, we recorded the 
electrical capacitance at 200 kHz and the electrical 
conductance at 100 Hz of the atmospheric ice between 
probe electrodes (Fig. 2) as a function of time both during 
and after ice growth, using a low-frequency impedance 
analyzer (HP4192A). The maximum sampling rate in our 
system was 0.5 Hz. We used copper-foil grid electrodes 
with a width and thickness of 0.68 mm and 0.075 mm 
respectively, spaced by 0.59 mm. The measurements were 
performed at temperatures ranging from -5 to -30 ®C, air 
velocities ranging from 15 to 115 km/h (or from 4.2 to 32 
m/s), and varying level of water content in air ranging from 
0.1 to 10 g/m3. 



(a) 



Time (sec) 


(b) 

Figure 3. Typical examples of ice-sample capacitance, C, 
and conductance, G, vs. time curves. The computed 
relation between Q and and that between and as 
calculated for solid ice are also shown for comparison, (a) 
Temperature = -10 ®C, air speed =115 km/h, water content 
in air = 1.67 g/m^, and ice growth rate = 13.25x10’^ mm/sec. 
(b) Temperature = -30 ®C, air speed =115 km/h, water 
content in air = 0.25 g/m^ and ice growth rate = 6.39x10"^ 
mm/sec. 


3. RESULTS AND DISCUSSIONS 

From our findings, figure 3 shows the typical time- 
variations of capacitance and conductance of the ice 
sample. These results not only show that the electrical 
properties of the simulated atmospheric ice changed 
significantly when the water spray was turned on and off, 
but also show a complex time dependence related to the 
experimental conditions. At a high rate of water flow (high 
air velocity times water content in air. Fig. 3a) and high 
temperature, the sample’s capacitance, C, monotonically 
decreased from a value typical for water to one more 
characteristic of solid ice; while at the same time, the 
sample’s conductance, G, monotonically increased 
(following a sharp rise upon turning on the water spray). In 
contrast, at a low rate of water flow and low temperature 
(Fig. 3b), the capacitance, C, gradually increased at a rate 
closely corresponding to the increasing solid ice thickness. 
During ice growth, however, the sample conductance, G, 
quickly decreased to a level characteristic of solid ice. 

Peaks seen on C- and G-curves in Fig. 3a were typical 
for tests with high rate of water flow after the water spray 
was shut down. One possible explanation of those peaks is 
that unfrozen water re-enters ice/grid electrodes interface 
when liquid inclusions in ice ultimately freeze, perhaps, 
due to ice fracture. 

Since the measured time-variation of the ice sample’s 
capacitance and conductance reflect not only LWCi but 
also the sample thickness (L^), we separate these effects 
using the relationship between C and and that between G 
and for solid ice. We determined these relationships 
using the electrostatic model and the conductive media DC 
model from a FEMLAB software package (Comsol Inc.). 
This software can solve the electrostatic boundary problem 
and the steady electric current problem by a finite element 
method. Figure 4 shows a schematic view of our model’s 
geometry. We performed a simulation of C and G for one 
unit cell and calculated the total capacitance and 
conductance of ice sample on the whole surface of the grid 
electrodes. Insulating boundary conditions are applied to 
all boundaries except the surface of the electrodes. The 
dielectric constant and the electric conductivity of ice are 
set as the values measured in our experiments for each test. 
Figure 5 shows an example of such simulation for the 
electric field in the ice using the electrostatic model. 

Figure 6 shows the C-L^ curve and the G-L^ curve 
obtained by the simulation. As seen in Fig. 5, the electric 
field penetrates in the ice only over a distance comparable 
with the inter-electrode space (« 0.6mm). Because of this, 
C and G do not change much after ice thickness exceeds 
that value, see Fig. 6. 
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Figure 4. Schematic geometry of simulation model. 



Figure 5. An example of simulation of the electric field in 
the ice. Ice thickness was 1.5 mm. Arrows show the electric 
field in the ice. The dielectric constant was set at 2.7. 



(a) 


- 100 

- 80 

- 60 

- 40 

- 20 

- 0 

I [ I I I I 

-100 0 100 200 300 400 500 

Time (sec) 



O 

_P 



Sample Thickness (mm) 


Figure 6. Relationship between the capacitance and the 
sample thickness and that between the conductance and the 
sample thickness obtained by a computer simulation for 
solid ice. The dielectric constant for ice was set at 2.7. The 
electric conductivity of ice was set at 1x10’^ S/m. 


(b) 

Figure 7. An example of the time-variations of C/C, and 
C/C, in growing atmospheric ice. (a) Temperature = -10®C, 
air speed =115 km/h, water content in air = 1.67 g/m^, and 
ice growth rate = 13.25x10'^ mm/sec. (b) Temperature = - 
30 ®C, air speed =115 km/h, water content in air = 0.25 
g/m^ and ice growth rate = 6.39x10'^ mm/sec. 

In Figs. 3a&b, the sample thickness effect on ice’s 
capacitance and conductance obtained by the simulation 
are plotted with the raw data for comparison, where the 
simulation results are plotted as a function of time 
assuming the ice growth rate is constant. The data 
normalized to the sample thickness effect can be expressed 
in terms of the ratio of C/C, and C/C„ where C and C are 
the sample’s capacitance and conductance respectively and 
C, and G, are the sample’s capacitance and conductance for 
liquid-water-fi*ee ice respectively (obtained by a 
simulation). In this way, the data indicate that LWCi 
approaches zero as these ratios approach one. Figure 7 
shows the ratios of C/C, and G/G, as a function of time. It’s 
important to remember that what is seen in Fig. 7 is 
essentially C and G of 0.6-mm thick pre-electrode layer of 
ice. 
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when it consists of all water and all ice, respectively. The 
total capacitance of the sample (Gtotai) for each connection 
model can be expressed in the same way as by 
replacing Q and with and G^, respectively. 


(a) parallel (b) series 


If we use the following notations: 


Y 1- Y 

HI wr X’ 

(c) combination 1 


Y l-Y 



1-Y 

X’ 


(d) combination 2 


Figure 8. Connection models between water and ice parts 
in the atmospheric ice. 


Rc=cjq 

Rg = GjG: 


(5) 


LWCi can be solved in terms of and Rg 

for each connection model as: 

(a) Parallel connection model 


To determine LWCi using results on ice capacitance 
and ice conductance, we have to assume some specific type 
of connections between frozen and imfrozen parts of the ice 
sample. Figure 8 shows four possible types of such 
connections: 


LWCi = X = 


^total ^ 

C,- y 


(^c-1) 


( 6 ) 


(a) parallel connection 

(b) series connection 

(c) combination #1 of parallel and series connections 

(d) combination #2 of parallel and series connections 


LWCi = X = 


G ^ 

^ total _2 

y 


(7) 


The total capacitance of the sample {C^otai) foi* ®^ch 
connection model can be expressed as: 


(a) Parallel connection model 
C^^,=XC^+{\-X)C^, 

(b) Series connection model 


( 1 ) 
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- + - 




CJY C,/(l-7), 
(c) Combination model #1 


( 2 ) 
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cjY^ q/(i-r). 


(d) Combination model #2 


C = 

'-'total 
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, (3) 
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C^/F (X'C, + (1-X')C,)/(1-F) 


(b) Series connection model 


C, 


total 


LWCi = Y = 


Ci 


R. 
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LWCi = Y = k 
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( 8 ) 


(9) 


(c) Combination model #1 

LWCi = X'+{\-X')T, (10) 

X’ and Y’ can be solved by the following equations with 
given Gfota/Gp Gf(,fa/Gi9 Rc 9 ^g- 


total 


RX\HRc-^)X'-iRc-^)X'T) 


( 11 ) 


(4) 

where X, 7, X' and 7' are fractions of water parts as shown 
in Fig. 8 and and Q are the capacitance of the sample 


R,-iR,-W' 
























































































(d) Combination model #2 


( 12 ) 

LWCi =(l-7')X'+7' 

X’ and Y’ can be solved by the following equations with 
given C^JCi, R„ and Rg. 

c,„., RA^ + (K-\)X') 

C, «,-(«,-i)r+(fi,-i)x'r 

> 

(14) 

g.., i?,(l + (i?,-l)X’) 

G, R^-iR^-X)T^iR^-\)TT 


The ratios of both and are given by the 

data of aCi and G/G^, The ratios of both CJC^ and GJG, 
are constants and can be determined experimentally. Here, 
Rc and Rg were set as 37 and 400, respectively. They were 
measured at 0 for water and -10 ®C for ice, respectively. 


Systems of equations (11)-(12) and (14)-(15) are 
quadratic for X' and Y' and, thus, return pairs of solutions. 
To choose a right pair we applied additional criteria: 
0 < X’< 1 and 0 < 7’< 1, and a requirement of minimum 
LWCi. 

Figures 9a&b show examples of the time-variations 
of LWCi as calculated with assumptions of the above four 
models. The results show that LWCi could vary depending 
on the connection model chosen. The combination models 
#1 and #2 have the similar results of LWCi. The reason is 
that the obtained X' and Y' showed the conditions of either 
X'«Y* or X*»Y' so that a term of X'Y' become smaller 
than terms of X' and Y\ As the results, the difference 
between two combination models diminished. 

The results also show that LWCi as calculated for the 
series model with use of the conductance data is 
unrealistically large and is not consistent with LWCi as 
calculated for the same series model with use of the 
capacitance data. Thus, we conclude that the series 
connection model is not suitable. 



Time (sec) 

(a) 



(b) 

Figure 9. Time-variations of LWCi as calculated in four 
models, (a) Temperature = -10 ®C, air speed =115 km/h, 
water content in air = 1.67 g/m^ and ice growth rate = 
13.25x10'^ mm/sec. (b) Temperature = -30 ®C, air speed = 
115 km/h, water content in air = 0.25 g/m^, and ice growth 
rate = 6.39x10'^ mm/sec. 



Figure 10. Thin section of atmospheric ice at high water 
rate conditions (75 km/h and 0.75g/m^) and -lO'^C. Air 
flow was parallel to elongating column grains. 











































We found a similar inconsistency of LWCi calculated 
from C- and G-data for the parallel-connection model, see 
Fig. 9. But either of two combination-connection models 
returned consistent and reasonable values of LWCi. For 
that reason, we give that model our preference. Moreover, 
examination of ice microstructure using thin sections cut of 
grown ice samples also support the assumption of mixed 
mode in which unfrozen water channels (grain boundaries) 
oriented in the electrode vicinity, see Fig. 10. Figure 10 
shows an example of thin sections of atmospheric ice at the 
high water rate (75km/h and 0.75g/m3) and -10 
obtained by our previous test. Small grains were formed at 
the bottom of the ice sample and large columnar grains 
were formed on the small grains. There are also isolated 
round grains between columnar grains. This picture 
indicated that we dealt with some sort of combination of 
series and parallel connections of liquid and solid regions 
in the ice samples. Therefore, we may conclude that a 
combination of series and parallel connections of liquid 
and solid regions is most appropriate MODEL. 

4. CONCLUSION AND FUTURE STUDY 

In order to measure liquid water content in growing 
atmospheric ice (LWCi) in-situ, we relied on the known 
differences in electrical properties between water and ice 
and recorded the time-variation of the capacitance and 
conductance of atmospheric ice under a wide range of 
experimental conditions. We developed a data analysis 
method of LWCi for the time-variation of the samples’ C 
and G, and then normalized the raw data for the sample 
thickness effect. LWCi could then be determined from the 
time-variation data for these electrical properties using a 
combination series and parallel connection model. 

In our analysis method, we ignored the dependence 
on temperature of the dielectric constant and conductivity 
(i.e., Rc and Rg). However, the time-variations of C and G 
also reflect the sample’s temperature changes during icing. 
Therefore, future study will involve taking this temperature 
dependence into account in order to further refine the data 
analysis method of LWCi. We will also perform similar 
experiments on grids of electrodes of various dimensions to 
investigate how LWCi varies with the distance from the 
interface. 
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